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Abstract 

The structure of the interstellar medium in the central kpc region of a galaxy with a weak bar-like 
potential is investigated taking into account realistic cooling and heating processes and the self-gravity 
of the gas. Using high resolution hydrodynamical simulations, it is revealed that the resonant structures 
(e.g. smooth spiral shocks and a nuclear ring) are very different from those seen in past numerical models 
where simple models of the ISM, i.e. non-self-gravitating, isothermal gas were assumed. We find that the 
pc-scale filaments and clumps form large scale spirals, which resemble those seen in real galaxies. The 
fine structures are different between the arms and in the nuclear region. The next generation millimeter 
interferometer (ALMA) may reveal the fine structures of the cold gas in nearby galaxies. 

We also find a large scale anisotropy in the gas temperature, which is caused due to non-circular velocity 
field of the gas. The damped orbit model based on the epicyclic approximation explains the distribution 
of the hot (> 10 4 K) and cold (< 100 K) gases appearing alternately around the galactic center. Because 
of the temperature anisotropy, cold gases observed by molecular lines do not necessarily represent the real 
gas distribution in galaxies. 

Position- Velocity diagrams depend strongly on the viewing angles. As a result, the rotational velocity 
inferred from the PV maps could be two times larger or smaller than the true circular velocity. 
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1. INTRODUCTION 

Gaseous response to the bi-symmctric gravitational 
potential, such as a stellar bar in galaxies, has 
been well studied numerically since early 70's (Miller, 
Prendergast, & Quirk 1970; Sorensen, Matsuda, & 
Fujimoto 1976; Matsuda & Nelson 1977; Huntley, Sanders, 
& Roberts 1978; Matsuda & Isaka 1980; van Albada 
& Roberts 1981; Schwarz 1984; Matsuda et al. 1987; 
Athanassoula 1992; Wada & Habe 1992; Friedli & 
Benz 1993; Wada & Habe 1995; Piner, Stone, & 
Teuben 1995; Heller & Shlosman 1994; Englmaier & 
Gerhard 1997; Fukuda, Wada & Habe 1998; Fux 1999; 
Patsis & Athanassoula 2000; Englmaier & Shlosman 2000; 
Fukuda, Habe, & Wada 2000). A rotating bar- like poten- 
tial causes resonances in gas orbits, and spiral-like den- 
sity enhancement or shocks are formed in the gas disk. 
However, almost all simulations and analyses in the past 
three decades assumed a simple model of the interstel- 
lar medium (ISM) with the isothermal equation of state, 
where the gaseous temperature is assumed a ~ 10 4 K. 
In some models, the energy equation is solved, but 10 4 
K cutoff is introduced (Friedli & Benz 1993). In an- 
other major numerical representation of the ISM, i.e. the 
"sticky particle" method or cloud-fluid model (Combes & 
Gerin 1985; Fukunaga & Tosa 1991), the ISM particles 



loose their energy via inelastic collisions, and the veloc- 
ity dispersion of the particles are kept roughly constant. 
The isothermal equation of state or constant velocity dis- 
persion would be a relevant approximation to investigate 
the global gas dynamics in galaxies. It is obvious, how- 
ever, that the interstellar matter is not an isothermal and 
uniform media (McKee & Ostriker 1977), the inhomoge- 
neous nature of the ISM should be taken into account to 
investigate gas dynamics in a local scale or on a galactic 
central region. In the case that we have no observational 
information about the ISM with fine enough spatial reso- 
lution (e.g. parsec scale) in external galaxies, the simple 
treatment of the ISM is enough to compare the models 
with observations. Hopefully we will have extremely high 
quality observations on the ISM in the next decade. High 
resolution observations with the next generation millime- 
ter/submillimcter interferometer (e.g. ALMA) will reveal 
the molecular gas structure with a hundred times finer 
spatial resolution. Therefore more realistic treatment of 
the ISM in simulations on a galactic scale is necessary to 
understand the ISM through comparison with the obser- 
vations. 

Recently Wada & Norman (1999, 2001) have investi- 
gated dynamics and structure of a self-gravitating gas disk 
in the central kpc region of a galaxy, taking into account 
realistic cooling (10 K < T g < 10 8 K) and heating. Using 
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high accuracy Euler mesh code, they showed that the glob- 
ally stable, inhomogcnoeus interstellar matter is formed as 
natural consequences of gravitational and thermal insta- 
bilities. The surface density ranges from 1CP 2 to 10 6 M Q 
pc~ 2 , and the dense gas formes clumps or filaments where 
the temperature is T g <~ 10 — 100 K. In the diffuse gas, on 
the other hand, the gaseous temperature is T g ~ 10 4 — 10 5 
K as a result of shock heating or T g ~ 10 6 — 10 7 K due to 
supernova explosions. 

In this paper, we invesigate the gas dynamics in a weak 
barred potential using the same numerical method pre- 
sented by Wada & Norman (2001) (hereafter WN01). 
Energy feedback from stars and supernovae, which are 
implimented in WN01, are not considered here for sim- 
plicity (see Wada, Spaans, & Kim 2001 for effects of the 
energy feedback to the ISM). 

This paper is organized as follows. In §2 we give a 
brief description of our numerical code and models. In 
§3.1 numerical results of the global structure of the self- 
gravitating multi-phase ISM in a weak barred potential 
are discussed. Effects of the advective heating and cool- 
ing are analytically studied in §3.2. In §4, observational 
implications derived from the position velocity diagrams 
and rotation curves are discussed. Finally we summarize 
the results in §5. 

2. NUMERICAL METHODS AND MODELS 

In order to model the global evolution of the multi- 
phase ISM, a hydrodynamical code must deal with a 
wide dynamic range (e.g., scale: 1 pc ~ a few kpc, den- 
sity: 10~ 3 ~ 10 6 cm' 3 , and temperature: 10 ~ 10 8 K). 
Moreover, shocks must be resolved correctly, because typ- 
ical Mach numbers range from 10 to 100 in the ISM of 
galaxies, and also SNe and stellar winds cause strong 
shocks with Mach numbers of several hundreds. We also 
need high spatial resolution for the hot, diffuse gas as well 
as for the cold, dense gas. In a rotating non-axisymmetric 
potential, the dynamical resonances play an important 
role for producing spiral shocks, density waves, and ring- 
like strucutres. We are interested in how the local inho- 
mogeneous structure of the ISM affects formation of these 
global resonant structures. 

We use hydrodynamical simulations for a single fluid 
taking into account the self-gravity of the gas, radiative 
cooling, and heating processes. The multi-phase feature of 
the gas, that is, the multi-temperature and multi-density 
structure, should be achieved as a result of the non-linear 
evolution of a single fluid, if we take into account realistic 
cooling and heating processes, and the self-gravity of the 
gas with relevant spatial resolution. We investigate the 
effects of the global dynamics of the gas in a bar poten- 
tial, therefore it requires that the whole rotating disk is 
simulated. The hydrodynamical part of the simulations 
is solved by a high-accuracy explicit Eulerian code with 
non-adaptive numerical grids. The effect of the radiative 
cooling is implemented using an implicit method with a 
cooling function. The numerical code is tested for various 
standard problems (see details WN01). 



We solve the following equations numerically on a ro- 
tating frame of a bar potential in two dimensions. 
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where, p,p,v are density, pressure, and velocity of the 
gas, the specific total energy E = \v\ 2 /2 +p/ (7 — l)p with 
7 = 5/3, and fij, = is the pattern speed of the bar 
potential. 

The potential of the analytical models is 



$ ext (R,cj>)=$o(R) + $ 1 (R,c 



(5) 



where $0 an d $1 are axisymmetric and non-axisymmetric 
potentials, respectively. We use the 'Toomre disk' poten- 
tial for the axisymmetric part. The potential has two pa- 
rameters a (scale parameter called 'core radius') and w max 
(the maximum rotation speed), and expressed as 

„2 



MR) = 



a(R 2 + a 2 ) 1/2 



(6) 



where c is given as c = v max (27/4) 1 ^ 4 a. The non- 
axisymmetric part of the potential is assumed to be in 
the form 



$i(i?, =e(R)$ cos2(<j)-n b t) , 
where e(R) is given as 



e(R) = £0 



aR 2 



(R 2 + a 2 ) 3 / 2 ' 



(7) 



(8) 



The parameter £ represents the strength of the bar. This 
formalism is taken from Sanders (1977). The parameter £0 
is at full strength from the beginning of simulations. Since 
our bar models are weak (e = 0.05, i.e. $i/$ ( a ) = 0.02 
) , initial strong shocks and subsequent accretion of gas to 
the galactic center do not occur. 

The equations (1), (2), and (3) are solved by second- 
order AUSM (Advection Upstream Splitting Method) 
based on Liou & Steffen (1993). We achieve third-order 
spatial accuracy with the Monotone Upstream-Centered 
Schemes for Conservation Laws (MUSCL) approach (van 
Leer 1977). 

Time integration is achieved by the second-order leap- 
frog method. Each time step is determined by the Courant 
condition. We adopt implicit time integration for the cool- 
ing term. 

Various two-dimensional test problems, such as the dou- 
ble Mach reflection of a strong shock, a point explosion in 
an adiabatic uniform medium, an isothermal gas flow in 
a weak barred potential, and evolution of spiral gravita- 
tional waves in a differentially rotating disk are described 
in Wada & Norman (2001). We find that AUSM with 
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MUSCL is as powerful a scheme for astrophysical prob- 
lems as PPM (Woodward & Colella 1984), but the algo- 
rithms and coding are much simpler than PPM. 

We use 2048 2 Cartesian grid points for calculating gas 
dynamics in a 2 kpcx 2 kpc region which includes the 
galactic center. Therefore the spatial resolution is 0.98 pc. 
The gravitational potential of the gas is calculated using 
the convolution method with the Fast Fourier Transform 
(Hockncy & Eastwood 1981). In this method, 4096 2 grid 
cells and a periodic Green's function arc used to calculated 
the self-gravity of the isolated matter in the 2048 2 grid 
points. 

We also assume a cooling function A(T S ) (10 < T g < 
10 8 K) and the photoelectric heating rate Tyy We use 
the cooling curves of Spaans & Norman (1997) (Fig.l in 
WN01). The cooling processes taken into account are, 
(1) recombination of H, He, C, O, N, Si and Fe, (2) col- 
lisional excitation of HI, CI-IV and OI-IV, (3) hydrogen 
and helium bremsstrahlung, (4) vibrational and rotational 
excitation of H2 and (5) atomic and molecular cooling due 
to fine-structure emission of C, C + and O, and rotational 
line emission of CO and H 2 . We assume solar metalicity 
in this paper. 

We assume a uniform UV radiation field and photo- 
electric heating of grains and PAHs. The heating rate 
Tuv is the same as in Gerritsen & Icke (1997): Tuv = 
1.0 x 10~ 24 eGo ergs s _1 , where the heating efficiency e is 
assumed to be 0.05 and Go is the incident FUV field nor- 
malized to the local interstellar value. 

The initial condition is an axisymmctric, uniform and 
rotationally supported disk, whose radius is 1 kpc, with 
the total gas mass 5 x 10 7 M Q which is about 2 % of the 
total mass determined by the background gravitational 
potential. The initial rotational velocity is chosen in order 
to balance the centrifugal force caused by the external 
potential <J> ext and $ sg . 

Random density fluctuations are added to the initial 
disk, with an amplitude less than 1 % of the unperturbed 
density and temperature. The initial temperature is set 
to 10 4 K over the whole region. In ghost zones at bound- 
aries, the physical values remain at their initial values 
during the calculations. From test runs, we found that 
these boundary conditions are much better than so called 
'outflow' boundaries, because the latter cause strong un- 
physical reflection of waves at the boundaries. 

3. NUMERICA1 RESULTS 

3.1. Non-isothermal, self- gravitating gas in a bar poten- 
tial 

For comparison with the realistic model, we first show 
the density and temperature distribution of a model with- 
out taking into account self-gravity of the gas in Fig. 1, 
which is a snapshot at t = 20 Myr (about 2.4 rotational 
period at R = 0.2 kpc). The surface density map (Fig. 
1(a)) shows a typical gas response to a weak bar poten- 
tial with two inner Lindblad resonances (IILR and OILR). 
The two inner leading spirals form an oval structure. The 
outer two-arm trailing spiral is a consequence of the OILR 



(Wada 1994). The temperature map (Fig. 1 (b)) shows a 
quadruple distribution. The gas is heated up to 10 4 ~ 10 5 
K at upstream sides of the leading and the trailing spi- 
rals. At downstream regimes, the gas temperature is as 
low as the lower limit of the temperature (10 K) to the 
contrary. This extreme anisotropy is because the advec- 
tive cooling/heating that dominates the heat balance in 
the supersonic non-axisymmctric flow. Sec details in §3.2. 

If we take into account self-gravity of the gas, how- 
ever, the density and temperature distributions in a quasi- 
stable state have very different fine structures from Fig. 1 . 
Figure 2 shows density and temperature maps of the self- 
gravitating model at t = 61 Myr. As seen in the density 
map (Fig. 2(a)), the clumpy and filamentary dense gases 
form weak spiral density enhancements, whose location 
roughly coincides with those in the non-self-gravitating 
model. The cold gases form complexes of clouds in the 
spirals that look similar to the giant molecular clouds in 
spiral galaxies. 

The complicated network of the filaments and clumps 
on a local scale is the same as that in axisymmetric mod- 
els in WN01. A number of processes are involved in the 
formation of the filamentary structure: 1) tidal and col- 
lisional interactions between dense regions formed due to 
the thermal and gravitational instabilities, 2) differential 
rotation, and 3) shear motion due to local turbulent mo- 
tion, 

The closeup of the central 500 pc region reaveals the 
fine structure of the gas (Fig. 3(a)). The dense clumps 
in the nuclear region are typically 20-50 pc across, and 
their masses ~ 10 5 Mq. The filamentary structures, which 
are nearly parallel to the bar major axis, are prominent. 
They are formed due to non-axisymmetric motion and 
shear in this region. The density range is greater than 
four orders of magnitude. On other hand, the closeup of 
the arm region (Fig. 3(b)) shows that its fine structure is 
very different from that in the nuclear region. The size of 
the dense clumps in the arms are smaller than that in the 
central region, because the velocity dispersion is smaller 
in the arm regions. This result suggests that the star 
formation processes in the arm and in the central regions 
would be different. 

Comparing the realistic ISM model (Fig. 2) with the 
non-selfgravitating model (Fig. 1), the continuous galactic 
shocks seen in the non-selfgravitating model, is no longer 
seen in the realistic model. Spiral structures are actually 
ensembles of dense clouds and filaments. The continuous 
central ring seen in Fig. 1 is not also reproduced in the 
self-gravitating model, where the central density structure 
is very complicated as seen in Fig. 3 (a). These results 
suggest that self-gravity of the gas should not be negligi- 
ble even if the gas mass fraction to the dynamical mass 
is small (2 % in the present case), in order to explore the 
gas dynamics in the central kpc region of galaxies. The 
isothermal approximation for the ISM, which was often 
used in past numerical simulations, would be relevant for 
a galactic scale gas dynamics, but more realistic treatment 
of the cooling processes should be necessary for the ISM 
on a smaller scale. We should take into account small 
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scale structure of the ISM especially for the cases that 
star formation and its energy feedback to the ISM arc im- 
portant. If one uses the isothermal gas model for the ISM 
with star formation, it is hard to avoid to introduce some 
phcnomcnological models for the star formation and en- 
ergy feedback processes. Our numerical approach, which 
is based on the basic equations that represent the phe- 
nomena, is less ambiguous in terms of these points. 

3.2. Effect of the advective cooling in a non-uniform 
stream motion in a bar potential 

We mentioned in the previous section that temperature 
of the gas in a bar potential changes quadruply on a global 
scale, i.e. cold (T g < 100 K) and hot (T g > 10 4 K) regimes 
appear alternately around the galactic center (Fig. 1(b) 
and 2(b)). This is because the advective cooling/heating 
dominates the heat balance in this system. This can be ex- 
plained analytically using the "damped gas orbit model" 
under the epicycle approximation. 

Wada (1994) provided an analytical representation 
which represents the behaviour of a non-self-gravitating 
gas in a rotating bar-like potential. Its outline is as fol- 
lows. The gas dynamics can be understood in terms of 
closed elliptical orbits which are solutions to the equations 
of motion of a forced-damping oscillator driven by a peri- 
odic external force. Gaseous elliptical orbits always incline 
to the bar potential in a leading sense inside of the corota- 
tion, and in a following sense outside of the corotation. If 
there are two ILRs, the gas orbits at the ILRs are oriented 
by 45° with respect to the bar potential. The damped or- 
bits gradually change their orientation with their radius. 
Since the direction of the gradual rotation near to the first 
ILR and that near to the second ILR are in an opposite 
sense, leading or trailing spiral-like enhancements appear 
around the first and the second ILR, respectively. This 
gradual rotation is explained in terms of the phase-shift of 
a damping oscillator. The analytical model well describes 
the gas behaviour seen in hydrodynamical simulations. 

Under the epicycle approximation, the simplest way to 
represent the effect of the gaseous nature, such as vis- 
cosity, on closed orbits is to introduce a simple damping 
term, which is proportional to the velocity of the radial 
oscillation 2A R\ into the equation of motion. 



Ri + 2\R 1 + k^Ri = / cos2(Q - ^ h )t, 



where 



dR 



2ft$v 



R(n-n h ) 



(9) 



(10) 



Ro 



and A is the damping rate. This equation describes a 
radial damped oscillation at around Ro, which is driven 
by an external periodic force at a frequency of 2(f2o — ^b)- 
The general solution to this equation is 

= Ae~ xt cos(ujt + a) + Bcos[2 (fi -fl b )t + S$l) 



where A and a are arbitrary constants, and u> = a/kq — A 2 . 
The first term of the right-hand side of (11) is a general so- 
lution of a harmonic oscillator having a natural frequency 



of uj with weak friction. Amplitude B and the phase-shift 
So = S(Ro), are 

fo 



^{k 2 - 4(ft - f! b ) 2 } 2 + 16A 2 (f> - fl b ) 2 

and 

Sq = arctan 



2FQ 



(12) 



(13) 



(14) 



^F 2 -l 

where 

F = 2(n Q -n h )/K Q , & = x/ Ko . 

The phase-shift So is always negative if F > 0, and So = 
— n/ 2 when F 2 = 1, that is, at the Lindblad resonances 
(flo = Ob ± Ko/2). The negative So means that the damp- 
ing oscillation is delayed from the periodic driving force. 
Since the first term on the right-hand side of equation (11) 
is a damping term which does not represent a closed orbit, 
we ignore it. Using (j>o(t) = (Qo — ^b) t, we obtain a closed 
orbit (Ri,(j>i), 

i?i(0o) = Scos(20 o + 5o) 

and 

fl B 



(15) 



01 (0o ) = - 



Ro(rio-tt b ) 



[sm(S ) -sin(20 o + 6o) 



MR) 



2R 2 (n -^b) 



■sin(20o). (16) 



We plot the orbits for the bar potential, 

$bar = (R 2 + & 2 r 1/2 [l + eobR 2 (R 2 + b 2 )- 3 ^ 2 



(17) 



with e = 0.05, b = 0.5, Vl b = 0.25, A = 0.1 in Fig. 4. 

The gas velocity field (vr,w</>) is obtained from (15) and 
(16): 



vr{Rq) 
f>(Ro,<fio) 



-B(Q - fib) sin(0 o + S ), 
[R + Bcos((j>o + So)} 
2O B 



Ro 



- cos(2<j(>o + S ) 



fl b {Ro) 



(18) 



R 2 (Q Q -flb) 



cos(2(^9) 



Now, the energy equation in the hydrodynamical simu- 
lations can be written as 



d P E dpH_ d P H 



pT vv - p 2 A(T g ) 



(20) 



dt dx dy 

where specific enthalpy H = E+p/p with E = \v\ 2 +p/ (7 — 
l)p. 

The advective cooling/heating rate is about one order of 
magnitude larger than the radiative cooling/heating rate 
in the supersonic gas disk that we are exploring. If the 
radiative cooling and heating is in an equilibrium (i.e., the 
right-hand side of eq. (20)), temperature change in a unit 
time due to the advective terms can be written as 

pm H (j- 1) 



AT ad (i?,< 



1500 



d\v\ 2 



d\v\ 2 d\v\ 2 



■ + 



dR 

d\v\ 2 



dR Rd(f> 



RB<f> 
KMyr" 1 



(21) 
(22) 
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where v 2 ^> c 2 and p ~ const is assumed for simplicity. 
Using (18) and (19), we plot - [d\v\ 2 /dR + d\v\ 2 / Rd(f\ 
in Fig. 5. AT ad is about ±10 3 - 10 4 K Myr" 1 . The 
quadruple anisotropy of the advective cooling and heating 
in the hydrodynamical simulations (Fig. 1 and 2) is well 
reproduced, and it suggests that the non-axisymmetric 
velocity field of the gas dominates the global temperature 
distribution in a weak barred potential. 

4. OBSERVATIONAL IMPLICATIONS 

From the numerical results shown in §3.1, we can in- 
fer how the molecular gas in nearby galaxies is observed 
with the submillimeter interferometers. It is crucial to 
treat realistically the cooling and heating processes of the 
ISM, especially when one tries to compare the numerical 
result with observations, such as molecular line observa- 
tions. As shown in §3.1, the temperatures of the ISM in 
a bar potential globally anisotropic. This is essential to 
interpret the observed molecular line intensity distribu- 
tion. To make detail comparison between the models and 
observations, one needs radiative transfer calculations for 
molecular lines (see Wada, Spaans, & Kim 2000 for the 
Large Magellanic Clouds). Alternatively, here we can use 
the density of the cold gas (T g < 100 K) component in- 
stead of the line intensity derived from detail line transfer 
calculations in order to infer the observed molecular gas 
distribution. 

Figure 6 (a) is the surface density distribution of the 
cold gas (< 100 K) with a 100 pc resolution. One may 
regard this map as a "CO integrated intensity map" of a 
galaxy in the Virgo cluster observed with the Nobeyama 
Millimeter Interferometer, for example. Two clumpy trail- 
ing spirals are clearly seen, and one might see a nuclear 
"pseudo ring" . There are dense regions at the base of the 
spirals and also at the nuclear region. Figure 6 (b) is the 
same as Fig. 6 (a), but with 10 pc resolution. It is clear 
that spirals are actually an ensemble of small clumps and 
filaments (see also Fig. 2(b)). Comparison of the two 'ob- 
servations' suggests that spatial resolution of the present- 
day interferometers are still too coarse to discuss small 
scale (e.g. 100 pc) structure of the cold gas in nearby 
galaxies. In other words, one should not be confused by 
the small scale inhomogeneity, which are comparable to 
the beam size, seen in observed intensity maps. 

The other important implication from our numerical 
results is that distribution of the cold gas derived from 
molecular line intensity (e.g. CO (1-0)) does not neces- 
sarily represent true mass distribution of the gas even on 
a galactic scale. The advective cooling and heating of 
the gas due to the non-circular velocity field in a bar po- 
tential dominate global distribution of the gaseous tem- 
perature. Suppose there is a uniform gas ring between 
the IILR and OILR. The temperature of the gas ring is 
roughly divided into two hot (T g > 10 4 K) and two cold 
(T g < 100 K) regions. Therefore we cannot observe gases 
in the hot regions with the CO (1-0) lines, for example, 
and non-axisymmetric distribution of the gas is appar- 
ently observed. A "hidden" hot gas component coexists 



with the nuclear "pseudo-ring" seen in Fig. 6 (a). 

In order to resemble the observed position velocity (PV) 
maps for the molecular gas, we have convolved the numer- 
ical result (Fig. 2) with different 'beam sizes' and velocity 
resolutions (Fig. 7 and 8). The PV maps obtained for 
three different viewing angles are shown. From compari- 
son between the three figures, one can find that the PV 
map strongly depends on the viewing angles (9 V ). Along 
the bar major axis (i.e. 8 V — 0°), about a two times larger 
line-of-sight velocity than the true circular rotational ve- 
locity (the dashed curve) is apparently observed. If we 
observe the same model along the minor axis of the bar 
(i.e. 6 V — 90°), the PV diagram represents nearly rigid ro- 
tation, and the apparent circular velocity is about a half of 
the real circular velocity at the core radius (R <~ 200 pc). 
Therefore if one observes the ISM in a edge-on galaxy that 
has a weak bar, and evaluate the dynamical mass from the 
rotation curve, then the mass could be four times larger 
or smaller than the real mass, depending on the view- 
ing angle. This viewing angle effect originating in the 
non-circular gas motion has been also known in previous 
isothermal simulations (Wada et al. 1994; Athanassoula 
& Bureau 1999). More quantitative analysis of the effect 
of non-circular motion on the rotation curves is discussed 
in Koda & Wada (2001). 

Figure 8 is the same as Fig. 7, but with low spatial 
and velocity resolution (100 pc and 5 km s _1 ) which re- 
sembles the observations with millimeter interferometers 
such as Nobeyama Milliliter Array. Although the surface 
density map does not reveal the fine structure of the cold 
gas, the low resolution PV map represents basically the 
same kinematical feature of the system. The most notable 
difference between the high and low resolution PV maps 
is that there are high 'intensity', steep velocity compo- 
nents (Av ~ 50 — 100 km s _1 ) in the high resolution map. 
These features are caused by the local non-circular mo- 
tion. Most prominent examples are seen in the central 
region (R < 100 pc) in the PV map (6 V = 0° and 45°) of 
Fig. 7. These features are originated from filamentary 
structures along the bar major axis seen in the density 
map. Such features cannot be resolved as discrete compo- 
nents in the low resolution PV map (Fig. 8). Therefore, 
in the low resolution PV maps, one might interpret the 
rotation curve steeply rises near the galactic center. One 
should note that even the bar potential assumed here is 
very weak ($i/<j>o < 0.02), the non-circular motion of the 
gas is large enough to cause the viewing angle effect on 
the PV diagrams as seen in Fig. 7 and 8. 

5. CONCLUSIONS 

The parsec-scale structure of the interstellar medium 
(ISM) in the central kpc region of a galaxy with a weak 
bar-like potential is discussed based on our two dimen- 
sional, high resolution hydrodynamical simulations, tak- 
ing into account realistic cooling and self-gravity of the 
gas. The results and conclusions are summarized as fol- 
lows. 
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• Taking into account realistic cooling/heating pro- 
cesses and self-gravity of the gas, the smooth spi- 
ral shocks which appear in the non-self-gravitating, 
isothermal gas model are not formed in the ISM. 

• The resonant structures due to the inner Lindblad 
resonances, such as spirals, have actually parsec- 
scale fine structures, which are filamentary and 
clumpy dense gases (§3.1; Fig. 2). The fine struc- 
tures are different between the arms and the nuclear 
region (§3.1; Fig. 3) 

• There is large scale anisotropy of the gaseous tem- 
perature, which is due to the non-circular velocity 
field of the gas (§3.1; Fig. 1 and 2). The damped 
orbit model based on the epicycle approximation ex- 
plains the quadratic distribution of the hot (> 10 4 
K) and cold (< 100 K) gases around the galactic 
center (Fig. 4 and 5). 

• Because of the temperature anisotropy, cold gases 
observed by molecular lines do not necessarily repre- 
sent all gas components even on a global scale (§3.2; 
Fig. 6). 

• Position- Velocity (PV) diagrams depend strongly on 
the viewing angles even if the bar component is very 
weak. The line-of sight velocity inferred from the 
PV maps along the major axis of the bar potential 
is about four times larger than that along the minor 
axis (§4; Fig. 7 and 8). 

• PV diagrams obtained using the same spatial and 
velocity resolutions in the present millimeter inter- 
ferometers globally show the same kinematics of the 
gas with ten times finer resolutions. However, the 
low resolution PV maps do not reveal the fine veloc- 
ity structure, and as a result, steeply rising rotation 
curves near the galactic center (R < 100 pc) can be 
apparently observed (§4; Fig. 7 and 8). 

• The spatial resolution of the present millimeter 
interferometers for external galaxies are not fine 
enough to discuss the gas dynamics on a hundred 
parsec scale, and we should wait for the next gener- 
ation radio interferometers, e.g. ALMA (Fig. 6). 

Most numerical simulations on the gas dynamics in a 
non-axisymmctric gravitational potential in the last three 
decades have assumed the simple model of the ISM. Under 
such approximation, the gas dynamics is well understood 
theoretically. However, we should investigate again the 
dynamics and structure of the ISM in galaxies using more 
realistic model based on the relevant basic equations. 

We are grateful to M. Spaans for providing us the cool- 
ing functions. Numerical computations were carried out 
on Fujitsu VPP5000 at the Astronomical Data Analysis 
Center of the National Astronomical Observatory, Japan. 
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Fig. 1. (a) Density distribution of a model without 
self-gravity of the gas. The major-axis of the weak bar po- 
tential of a pattern speed, Ut = 140 km s _1 kpc~ 2 is horizon- 
tally fixed. 512 X 512 grid points are used. Color represents 
log-scaled density between 1 and 1O 3 M0 pc -2 . Gas rotates 
counter clockwise. Inner Lindblad Resonances are rocated at 
R = 0.2 and 0.4 kpc. (b) Same as (a), but for temperature 
between 10 and 10 5 K. 



Fig. 4. Gas orbits under the influence of the two ILRs, de- 
rived from the damped orbit model (eqn. (15) and (16)). 



Fig. 2. Same as Fig. 1, but self-gravity of the gas is taking 

into account. 2048x2048 grid points are used. Fig. 5. Distribution of adiabatic cooling/heating rate, 

-[d\v\ 2 /dR + d\v\ 2 /Rd<j>\. See eqn. (22). Contours are 
drawn every 1 from -10 to 10. Solid/dotted lines represent 
heating/cooling regions. 

Fig. 3. (a) Close-up of the central 500 X 500 pc region in the 
same model of Fig. 2. The color scale represents log-scaled 
density between 0.1 and 1O 3 M0 pc~ 2 . (b) Arm-region 



Fig. 6. (a) Density distribution of the cold gas (T g < 100 K) 
in a weak bar potential convolved with a 100 pc "beam" size, 
(b) Same as left, but with a 10 pc beam size. Unit of scale is 
kpc. 



Fig. 7. Surface density distribution of the cold gas with 10 pc 
grid size and 1 km s _1 for velocity resolution. Position- velcity 
diagrams from three viewing angles (0, 45, and 90 degree 
from the bar major axis) are plotted. The dashed curve is 
a rotation curve derived from the axisymmetric gravitational 
potential cq. (6). 



Fig. 8. Same as Fig. 7, but with 100 pc and 5 km s 1 spatial 
and velocity resolution. 
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